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ABSTRACT 

We report on the first SCUBA detection of a Type 2 QSO at z=3.660 in the Chandra 
Deep Field South. This source is X-ray absorbed, shows only narrow emission lines in 
the optical spectrum and is detected in the sub-mm: it is the ideal candidate in an 
evolution scheme for AGN (e.g. Fabian (1999); Page et al. (2004)) of an early phase 
corresponding to the main growth of the host galaxy and formation of the central black 
hole. The overall photometry (from the radio to the X-ray energy band) of this source 
is well reproduced by the spectral energy distribution (SED) of NGC 6240, while it 
is incompatible with the spectrum of a Type 1 QSO (3C273) or a starburst galaxy 
(Arp 220). Its sub-mm (850 /zm) to X-ray (2 keV) spectral slope (asx) is close to the 
predicted value for a Compton-thick AGN in which only 1% of the nuclear emission 
emerges through scattering. Using the observed flux at 850 /jm we have derived a 
SFR=550-680 M /yr and an estimate of the dust mass, M dust = 4.2 x 10 s M Q . 
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1 INTRODUCTION 

The presence of massive black holes (MBHs) at the cen- 
ters of most galaxies appears by now firmly established. 
A number of correlations have been observed between 
the masses of MBHs and the properties of the g alacti c 
bulges hosting them (e.g . iKormendv fc Gebhardd 1.2000), 
iMerritt fc Ferraresd <l200ll) ). Since these correlations extend 
well beyond the direct dynamical influence of the MBH, it 
seems likely that a close link exists between the growth of 
the MBH and the formation of their host galactic bulges. 

If these two events are co-eval then we would expect 
that powerful black holes are hosted in the centers of the 
most massive galaxies. SCUBA galaxies show a strong red- 
shift clustering thus indicating that t hey are hosted by 
very high mass halos felain et alJ i200St) ). This is confirmed 
by CO molecul ar gas emissi o n lin e widths and possible 
rotation curves iFraver et alJ dl99gft : iGenzel et all J2003lk 
iNeri et all (120031) ). The combination of deep X-ray obser- 
vations and spectroscopic follow -up campaigns of SCUBA 
sources ^Alexander et all i2004f) ) led to the discovery that 
at least 38% of the SCUBA sources host an AGN, although 
in the majority of the cases the contribution of the nucleus 
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is not bolometrically significant (i.e., < 20%). 
A possible scenario would predict that the black hole grows 
within a SCUBA galaxy until t he obscu r ing g as is blown 
away revealing an AGN (e.g. iFabian 1 (1999)). A simi- 
lar evolution argument h as been proposed for ULIRGs 
iSanders fc Mirabell dl996tl ). which can be considered the 
closest local analogs to SCUBA galaxies. We have started 
to investigate this scenario using a unique sample of heav- 
ily X-ray absorbed (log(Ne) > 22 cm -2 ) and very luminous 
(log(Lx) > 44 erg s" 1 ) AGN, the so-called Type 2 QSOs, 
located in the Chandra Deep F ield-South (CDF-S) and se- 
lected from lSzokolv et all (120041) . In this letter we report on 
the first sub-mm detection of a high redshift (z=3.660) Type 
2 QSO and discuss the implications for its nature. We will 
discuss in a future paper the properties of the full sample 
of Type 2 QSOs once our observational programme is com- 
pleted. 

We assume a cosmology with parameters Qm =0.3, Qa = 
0.7 and H = 70 km s" 1 Mpc" 1 . 



2 SUB-MM OBSERVATION 

The Submilli metre Commo n User Bolometer Array 
fSCUBA. lHolland et alJ (|l9J)9j)) was employed in photome- 
try mode, with the wide 850:450 filter set, and a standard 60 
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Figure 1. Left: cutout in the z (F850LP) ,4 OS band of 11" across. 
The big dashed circle is the 5-arcsec radius from which the 850 
lira flux can be originated, while the smaller continuum circle in 
the center indicates the Chandra 3 a positional error. Right: a 
zoom of the central part (5" across) in the ISAAC K s band, the 
isointcnsity X-ray contours at 3, 11, 12, 13, 14 a above the local 
background are shown. 



arcsec chop in azimuth at 7.8Hz. The source was placed in 
the central bolometer (H7) and the median of the remain- 
ing bolometers was used for additional sky removal. Flux 
calibration was performed using Uranus. Telescope point- 
ing was checked frequently while sky opacity was monitored 
via regular skydips at 850 and 450/^m and continuosly via 
the JCMT Water Vapour Monitor and the CSO Tau Me- 
ter. T225 ranged from 0.05 to 0.13. The data were reduced 
using the SURF (SCUBA User Reduction Facility) software 
package. A linear interpolation between consecutive skydip 
values was used, and the nearest calibration source in time 
applied to account for the gain. After sky subtraction, the 
bolometer time-stream was clipped at the 3-cr level. Each 
of the individual datasets have been tested for consistency 
with one another using a Kolmogorov-Smirnov test, reject- 
ing anything below the 5% mark. 

We have a firm detection (S/N~ 4) for CDFS-263 1 
(RA= 3 : 32 : 18.86 DEC= -27 : 51 : 35.5, z= 3.660) : 
Sg50 = 4.8 ± 1.1 mjy. Since we already know the position 
of the X-ray source to sub-arcsec accuracy, the ~ 14 arc- 
sec FWHM SCUBA beam would imply that the sub-mm 
emission is either coming from this source or something else 
within a radius of ~ 5 arcsec. Inside this area we find two 
optical counterparts in a deep ACS z band image centered 
on CDFS-263 (see left panel of Fig. EJ. The X-ray contours 
clearly point to the counterpart named 'B' in Fig. Q (the 
continuum circ l e show s the 3cr Chandra positional error from 
iGiacconi et al] ]2002l) ) . The other possible sub-mm coun- 
terpart 'A' has a complex sub-arcsec structure, where all the 
different components have the same optical/near-IR colours 
therefore likely are at the same redshift. Since a good spec- 
troscopic redshift is not available for this object, we have 
derived a photometric redshift of z p hot = 1.6 ± 0.1 using 
the photometry in seven optical/near- IR bands (Tab le 3) 
and the public avail able code BPZ of iBenitezl i2000T) (see 
M ainieri et al] ^2004) for details on the procedure). Assum- 
ing this value, source 'A' has to have Lfir ~ 10 1 to shine 
in the sub-mm, and the probability to have by chance such a 
powerful source in the SCUBA error circle is low. The same 
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Figure 2. VLT FORS spectrum of CDFS-263 (z= 3.660) from 
Szokoly et al. (2004). 



Table 1. Optical emission line properties of CDFS-263. FWHM 
and EW are rest frame 



line 


flux 


FWHM 


EW 




[10~ 18 erg cm" 2 s _1 ] 


[km s _1 ] 


[A] 


Lya A1216 


2.96±0.04 


1000±10 


29.6±3.0 


N V A1240 


0.17±0.03 


1580±140 


2.4±0.6 


C IV A1548 


0.40±0.03 


1400±50 


5.7±0.9 



XID from lGiacconi et al] <200: 



area of the sky imaged in the Ks band with IS A AC /VLT 
shows only the 'B' counterpart (right panel of Fig. rein- 
forcing its chances to be the real sub-mm source. 



3 OPTICAL AND X-RAY SPECTROSCOPY 

To help assign a reliable classification of CDFS-263 as a 
Type 2 quasar we have determined its optical emission line 
and X-ray spectral properties. 

VLT FORS1 spectra of CDFS-263 were taken in 2000 
during the spectroscopic identifi cation programme of the 
CDF-S field JSzokolv et al] ]2004l) ). The object was observed 
with a multi-slit mask using 1.2 arcsec wide single slits and 
grism 1501 with a spectral dispersion of ~ 5.5 A/pixel, which 
corresponds to a spectral resolution of ~ 11 A. The seeing 
was about 0.5-0.6 arcsec. The total integrated exposure time 
of the spectrum (see Fig. |2J is about 3.3 hours. The mea- 
sured redshift is z = 3.660 ± 0.005. 

The optical spectrum of CDFS-263 shows a strong 
LycvA1216 emission line and the two fainter lines of N V 
A1240 and C IV A1548, whereas no further lines as e.g. 
Ly/3A1026 or SI IV A1397 are visible. Each emission line was 
fitted with a single Gau ssian profile a pplyin g the Levenberg- 
Marquardt algorithm dPress et al] ]l992h ). The four ad- 
justable parameters are the total line flux, the central wave- 
length, the FWHM in A and the flux of a local linear con- 
tinuum. From this set of parameters we have calculated the 
FWHM in km s _1 (corrected for instrumental resolution) 
and the rest frame EW in A. The line flux, the rest frame 
FWHM and the rest frame EW of the visible lines are given 
in Table 1. 
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Table 2. Best fit spectral model parameters for CDFS 263 
(9.3596E+05 sec) Galactic N H = 9.13 X 10 19 cm -2 . If we include 
the galactic N^j there is only a marginal change of the intrinsic 
Njj value. 



Chondra spectrum 



src_1M_263_g10.fits 



parameter 



PL + intrinsic absorption 
modcll model 2 



N H (10 22 


cm 2 ) 


Q +140 
u -0 


B8±? 5 


z 




3.660 (froozen) 


r 




34+ 2 -59 


1.90 (frozen) 


F p 15 




1 04+ - 66 
1 - u ^-0.52 


, 71 +1.37 
L - ,1 -3.42 


X 2 /dof 




1.19/4 


3.37/5 


fo. 5-10.0 


(absorbed) 


1.93E-15 


1.04E-15 


fo. 5-10.0 


(unabsorbed) 




2.00E-15 


Lo. 5-10.0 


(absorbed) 


0.19E-44 


0.37E-44 


Lo. 5-10.0 


(unabsorbed) 




7.64E-44 



Fit parameters are shown with 90% errors (Ax 2 = 2.7). The 
model and parameter definitions are: Model Components- PL: 
Power-law with a photon index of T and a 0.5-7 keV flux of F p i$ 
in units of 10 — 15 erg cm -2 s -1 using the XSPEC model pegpwrlw. 
Intrinsic absorption: Photo-electric absorption using Wiscon- 
sin cross-sections zwabs (Morrison & McCammon 1983), where 
N^f is the equivalent Hydrogen column in units of 10 22 atoms 
cm -2 and z is the redshift. 

The luminosities reported in the table are in the rest frame [0.5- 
10] keV energy band. 



N V A1240 is observed at the 4a level. The line fluxes 

are a factor of ~ 20 100 lower c ompared to tho s e of t he 

prototype Type 2 QSO CDFS-202 jNorman et alJ J2002l) 1. 
The emission line widths are narrow with a FWHM of < 
1500 km s _1 , whereas the mean line width of the permitted 
lines for X-ray selected Seyfe rt 1 and QSOl range s from 
~ 3000 to ~ 5000 km s" 1 (see iLehmann et all fcOOll) ). The 
narrow line widths and the observed 0.5-10 k eV luminosity 
of ~ 1 .6xl0 44 erg s _1 indicate a Type 2 QSO JSzokolv et alJ 
J2004h ). 

The analysis of the X-ray spectrum of CDFS-263 is es- 
sential to determine its amount of intrinsic X-ray absorp- 
tion and its unabsorbed X-ray luminosity. We have used the 
public Chandra ACIS -I data of the 1 Msec CDF-S survey 
iGiacconi et alJ i2002l) ). The total integrated exposure time 
is about 936 ksec for CDFS-263. The ACIS-I spectrum of the 
source was extracted from a circular region of ~ 5 arcsec in 
radius and the background spectrum was obtained from a 
source free region. The source spectrum contains 120 counts. 
The spectra were grouped to have at least 10 photons per 
bin. 

Due to the small number of source photons we can only 
get an estimate of the X-ray spectral properties. The spec- 
trum was analysed using XSPEC 11.3. A simple absorbed 
powerlaw model (1 in Table 2) gives a reasonable fit. The 
photon index of T = 0.34 is low compared to that of F ~ 1.9 
found for most Seyfert galaxies and quasars. This points to 
a large intrinsic absorption of CDFS-263. Therefore we have 
frozen the photon index to 1.9 (model 2). The results of 
both models are shown in Table 2. All the errors quoted we 
calculated at 90% confidence level for two interesting pa- 
rameters. As expected, the best fit model 2 shows a large 
Nh of ~ 10 24 cm -2 , which implies nearly Compton-thick 
absorption. If we correct for intrinsic absorption the 0.5-10 




channel energy (keV) 



Figure 3. The observed Chandra ACIS-I spectrum of CDFS 263 
and the best- fit absorbed powerlaw model (2). The spectrum was 
grupped to have at least 10 photons per bin. The bottom panel 
shows the ratios between the observed data points (crosses) and 
the best- fit model data points (continuum line). 



keV rest frame luminosity of CDFS-263 increases by a fac- 
tor of ~ 20 (see Table 2). The best fit model (2) and the 
X-ray data points are shown in Fig. |^1 Whereas the quality 
of the ACIS-I spectrum is low, the strong intrinsic absorp- 
tion of CDFS-263 (at least N H > 1 x 10 23 cm" 2 ) is fairly 
well constrained. 

The narrow optical emission lines (FWHM< 1500 km 
s _1 ), the large intrinsic absoprtion (log Nh > 22) and the 
large absorption corrected X-ray luminosity (log Lx > 44) 
in the 0.5-10 keV energy band favour the Type 2 quasar 
nature of CDFS-263. 



4 DISCUSSION 

4.1 Spectral energy distribution 

CDFS-263 has been observed in a wide range of energy 
bands. Photometry of extraordinary quality in the BVIz 
filters has been obtained with the ACS c amera on board 
of HS T as part of the GOODS survey iGiavalisco et al] 
(2004)) 2 . As part of the same survey, ESO carried out deep 
near-IR imaging with ISAAC/VLT in J, H and K s bands 3 . 
Finally, we have a radio upper limit from VLA observa- 
tions at 20 cm of F„ < 100 /iJy (K. Kellerman, private 
comunication) . The available photometry for CDFS-263 in 
summarized in Table 3. 

Accor ding to the radio-power criterion given bv lStocke et alJ 
(1992), we can classify this AGN as radio-quiet. Hence, a 
non-thermal (synchrot ron) contribution to the sub-mm flux 
is exp ected to be small ijCarleton et alJ l|l987f ): lrlughes et alJ 
jl993ft 1. 

We compare the photometry of CDFS-263 with the SED 
of several well-studied local sources. We choose 3C273 as a 
template of an unabsorbed quasar, NGC 6240 a well-studied 



2 These data are publicly 
|http : / / www . stsci . edu/science /goods/ 



These data are pub icly 

|http://www.eso .org/scicnce/go ods/| 



available at: 
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Table 3. Photometry of CDFS-263. 

Observed Rest-frame Observed vF„ 

Band Band AB Magnitude [W/m 2 ] 




Restframe Frequency [Hz] 



2-10 keV a 


27.9 keV 




2.6E-19 


0.5-2 keV a 


5.8 keV 




8.1E-20 


F435W 6 


927A 


27.7±0.4 


2.9E-19 


F606W b 


1277A 


25.67±0.05 


1.4E-18 


F775W b 


1654A 


25.38±0.08 


1.3E-18 


F850LP b 


1957A 


25.08±0.08 


1.6E-18 


J c 


2690A 


24.5±0.1 


1.4E-18 


H d 


3550A 


23.9±0.4 


1.8E-18 


K s c 


4660A 


22.46±0.04 


5.0E-18 


850/xm e 


1.6E3 GHz 




1.6E-17 


20 cm f 


6.9 GHz 




<1.5E-21 



a Chandra/ ACIS; b HST/ ACS; c VXT/ISAAC; d NTT/SOFl; 
e JCMT/SCVBA; f VLA 



Figure 4. Detections and upper limits for CDFS-263 plotted 
over the SED of NGC 6240 {solid), Arp 220 (dotted) and 3C273 
(dash- dotted). The dashed line shows where the photometry is not 
available and template is extrapolated. The CDFS-263 data have 
been shifted to the rest frame and the empty circles show the line 
subtracted magnitudes (see text). 



(U)LIRG which hosts a Compton-thick AGN, and Arp 
220 as a starburst galaxy where X-ray binaries are t he 
major source of its X-ray emission (Iwasaw a et al] J2001)). 
Due to the difference in redshift between our source and 
these templates, we have photometry for CDFS-263 in 
wavelength regions not covered for the local sources. In 
particular for CDFS-263, the V filter contains Lya while 
the i band contains the CIV line. We have subtracted 
the contribution of these lines in the observed fluxes and 
compared them with the templates SED in Fig. [I] The 
wide-band energy distribution of CDFS-263 is inconsistent 
with the spectrum of an unabsorbed QSO (dash-dotted line 
in Fig. : if we anchor 3C273 to the X-ray observations, we 
underpredict the fluxes at the other wavelengths. Using the 
SED of NGC 6240, we can reproduce in a consistent way 
the observed sub-mm, optical, near-IR and X-ray fluxes 
of our source (solid thick line in Fig. |1J. The observed 
flux in the B band is low compared with the expected 
one, this can be ascribed to absorption due to the Lya 
forest for a source at z = 3.660 (see also IStern et al] 
J2002I) ). Finally, the energy distribution of CDFS-263 does 
not resemble the SED of Arp 220 (dotted line in Fig. 0J: 
if we fix the template at the observed sub-mm flux we 
underestimate the X-ray flux by two orders of magnitudes. 
We note that Arp 220 appears to be an order of magnitude 
underluminous in soft X-rays, given its FIR luminosity, 
if compared w i th oth er starbursts galaxies (see Fig. 3 of 
llwasawa et al] i200 ll) ). but this is not enough to recon- 
cile the photometry of CDFS-263 with a starburst template. 



4.2 Sub-mm to X-ray spectral index 

To further compare the properties of our source with local 
templates and other SCUBA observations of X-ray sources, 
we derive the sub-mm (850 /im) to X-ray (2 keV) spectral 
slope (asx) of CDFS-263. Fig.|^|shows the expected values 
of asx, as a function of redshift, for a number of template 
SEDs. Namely, a quasar with no absorption in the X-ray or 
with different amounts of absorption (Nh = 10 23 , 10 24 and 



Compton-thick with 1% of the nuclear emission scattered); 
a mean starburst templ ate and Arp 220. Th ese templates 
are described in detail in lAlmaini et alJ l|2003F) . 
CDFS-263 yields a value of asx = 1-29 (star symbol in 
Fig. This value is clearly incompatible with an unab- 
sorbed AGN, as we already expected from the optical and 
X-ray spectral properties of the source. On the other hand, 
it is close to the spectral index of a Compton-thick AGN 
in which only 1% of the nuclear emission emerges through 
scattering. The X-ray spectral fit yields a column density of 
~ 10 24 cm -2 (see Table 2). One way to obtain this value of 
asx is a sub-mm flux due to starburst activity plus an ab- 
sorbed AGN in the center of the host galaxy that accounts 
for the X-ray emission, which is higher than that expected 
for a starburst gala x y alone (dashed l i ne in Fig. \5j (s ee also 
iFabian et all <200fJ) : lAlexander et all <2003ft : lAlmaini et al] 
(2003)). For comparison we report in Fig. [3 other X-ray se- 
lected sources with sub-mm data measurem ents from the lit- 
eratur e. Namely, the CDF-N sources from I Alexander et al] 
( 2003) (circles) the majority of which show asx values higher 
than CDFS-263, and compatible with the expected indices 
for a starburst galaxy, although all except one of them lack 
a precise redshift measurement. The X-ray absorbed (filled 
triangles) and X-ray una b sorbe d (e mpty triangl e ) R OSAT 
BLAGN from lPage et al] (1200 ll) and lPaee etaJ <|2004|) are 
sampling a different region in Fig. |H| with their asx val- 
ues being closer to pure QSO tracks in agreement with 
their nature as unabsorbed AGN in the optical band. Fi- 
nally, we show also the asx values of some z > 4 Type 1 
QSOs (square) with X-ray llKaspi et al] ((2000) ) and sub- 
mm IIMcMaho n et al] 11199 9) ) observations. The location of 
our source in the asx vs z plot is peculiar: it has a lower 
asx than any source from lAlexander et al] <|2003 | ) at high 
redshi ft and a higher value of the BLAGN from lPage et al] 
J200ll) fc lPaee et aH J2004I) and the Type 1 QSOs at z > 4. 
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given bv lMarconi et all i2004h estimate the bolometric lumi- 
nosity, Lboi = 2.4 x 10 46 erg s _1 . Assuming that Lboi =LEdd, 
where LEdd is the Eddington luminosity, we derive a lower 
limit on the black hole mass: M BH = 1.9 x 10 s M . 



5 CONCLUSIONS 

We have a firm (4cr) SCUBA detection of a Type 2 QSO at 
redshift z = 3.660 (CDFS-263). Its SED is consistent, from 
radio to X-ray wavelengths, with the shape of NGC 6240 
while it can not be reproduced with the SED of 3C273 or 
Arp 220. Futhermore, we derive an qsx = 1.29. We conclude 
that CDFS-263 can be powered by a Compton thick AGN 
(we found N H ~ 1 - 20 x 10 23 cm" 2 ). 

CDFS-263 is a strong candidate for an AGN in the ini- 
tia l phase (be f ore th e "X-ray absorbed phase") described 
bv lPaee et al.1 (2004) that corresponds to the main growth 
period of the host galaxy spheroid. It is a luminous X-ray 
source (L x > 10 44 erg s" 1 ) with high X -ray absorption 
(Nh > 10 23 cm" 2 ) and detected in the sub-mm, Ssso = 
4.8 ± 1.1 mJy. 



Figure 5. Sub-mm to X-ray spectral index (ogx) as a function 
of redshift. The dotted and dashed lines are the expected values 
of a sx for a set of SEDs (sec text). The values for CDFS-263 
(star), the sample from Alexander et al. (2003) (circles), Page et 
al. (2001) (filled triangles), Page et al. (2004) (empty triangle) 
and a sample of Type 1 QSOs at z > 4 (squares) are shown. 



4.3 Star Formation Rate, dust mass and black 
hole mass 

Under the hypothesis that the sub-mm flux is due to star for- 
mation activity, we estimate, using the measured flux at 850 
jUm, the dust-enshrouded Star Formation Rate (SFR) and 
the dust mass. We note that the temperature of the dust 
heated by the AGN would be too high (> 100 K) to emit in 
the sub-mm. The far-infrared luminosity (Lfir) provides a 
reasonable measure of the active formation of massive stars. 
Assuming a dust temperature o f 40K and an emissivity in- 
dex (3 — 1.2 (e.g. iDunne et all (I2000T0 we interpolate from 
the 850 /J,m luminosity and get an estimate of the 60 fim 
luminosity assuming an optically-thin greybody. The SFR 
is given by: 

SFR = elO- 10 L 60 /L aun (M Q /yr) 
where e is the fraction of optical/UV light being radiated in 
the far-infrared and depends on the IMF and the fraction 
of the optical/UV light that is absorbed and re-radiated in 
the FIR. For e = 0.6- -2 we get SFR = 550 - -68OM /yr 
(assuming Q\ = 0.7 and Ho = 70). 
The dust mass is given by: 

M duat = 1/(1 + z)Szs,oD 2 L /k d ( r «, st) B(vrest,T d ) 
where z is the redshift, Ssso is the observed flux density, 
Dl is the luminosity distance, k res t is the rest frequency 
absorption coeffi cient and B i s the Planck function. Using 
K re3t = 0.15 fsee lScott et alJ i2002t) ') we derive an M du3t = 
4.2 x 10 8 Mq. The dust mass estimate assumes optically thin 
thermal emission with no contribution from bremsstrahlung 
or synchrotron radiation. 

Finally, we compute the rest frame absorption corrected 
Lx[2 — 10 keV] = 4.5 x 10 44 erg s _1 and using the relation 
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